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Abstract 

Complex numbers appear in the Hilbert space formulation of quan- 
tum mechanics, but not in the formulation in phase space. Quantum 
symmetries are described by complex, unitary or antiunitary opera- 
tors defining ray representations in Hilbert space, whereas in phase 
space they are described by real, true representations. Equivalence 
of the formulations requires that the former representations can be 
obtained from the latter and vice versa. Examples are given. Equiva- 
lence of the two formulations also requires that complex superpositions 
of state vectors can be described in the phase space formulation, and 
it is shown that this leads to a nonlinear superposition principle for 
orthogonal, pure-state Wigner functions. It is concluded that the use 
of complex numbers in quantum mechanics can be regarded as a com- 
putational device to simplify calculations, as in all other applications 
of mathematics to physical phenomena. 
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1 Introduction 

One of the remarkable features of quantum mechanics as usually formulated 
is the fundamental role played there by the complex numbers. This ap- 
pears to distinguish quantum mechanics from other mathematical models of 
natural phenomena. Theoreticians use complex numbers routinely in most 
applications of mathematics, starting with the solution of simple quadratic 
equations. However, in all other models, the complex numbers are introduced 
only as a computational tool. Not only are the 'observables' of these models 
real but, invariably, so also are their defining equations. 

Feynman and Dirac were prominent among those who have stressed the 
crucial role of complex probability amplitudes in quantum mechanics. 

To quote Feynman: It has been found that all processes so far observed can 
be understood in terms of the following prescription: To every process there 
corresponds an amplitude (a complex number); with proper normalization the 
probability of the process is equal to the absolute square of the amplitude jl]. 

And Dirac: So if one asks what is the main feature of quantum mechanics, 
I feel inclined now to say that it is ... the existence of probability amplitudes 
which underly all atomic processes. Now a probability amplitude is related to 
experiment but only partially. The square of its modulus is something that 
we can observe. That is the probability which the experimental people get. 
But besides that there is a phase ... (which) is all important because it is the 
source of all interference phenomena but its physical meaning is obscure j2] . 

The essential place of complex numbers in quantum mechanics as usually 
formulated is evident in Schrodinger's time-dependent wave equation 

**2> = km . (i) 

where ip(t) is the state vector of a quantum system and H is its Hamiltonian 
operator. For spinless, nonrelativistic systems, H is a real operator in the 
usual coordinate realization of Hilbert space, and the radical departure that 
the complex equation (0) then represents, from the reality of the governing 
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equations of all earlier models of natural phenomena, was a departure that 
Schrodinger made uneasily [3]. 

There is a formulation of quantum mechanics jl], [3], [H] — I will refer to 
it as the phase-space formulation - that has become popular as an arena in 
which to investigate such matters as quantum chaos [7j, quantum tomogra- 
phy [S], 0, the relationship between quantum and classical mechanics [TU] . 
[11], and the nature of quantization [T2], [131 ■ What I would like to draw 
attention to here is that the phase-space formulation is purely real - there 
are no complex numbers to be seen in the defining equations. Despite this, 
it has been stated [12], [Mj that the phase-space formulation is equivalent 
to the more familiar formulation in terms of hermitian operators acting on a 
complex Hilbert space. In view of the remarks above, one may well ask how 
this can be so. 

In what follows, I shall discuss two aspects of this surprising state of 
affairs. The first concerns representations of symmetry groups and alge- 
bras. In the Hilbert space formulation, elements of a symmetry group or 
more generally, of any group of automorphisms of the set of quantum states, 
are represented by unitary or antiunitary operators that define ray repre- 
sentations of the group as a whole, in accordance with Wigner's celebrated 
theorem ^SJEIIE]- I n contrast, in the phase-space formulation, elements 
of a symmetry group are represented by real unitary operators that define a 
true representation. If the phase-space formulation is equivalent to the one 
in Hilbert space, then the complex ray representations must be recoverable 
from the real true ones, and vice versa. That this is indeed the case JH] 
is remarkable from the point of view of the mathematics and the physics 
involved. 

The second aspect concerns the fact that in the phase-space formulation, 
quantum states are described by real Wigner functions, whereas in Hilbert 
space, one has state vectors that can be superposed, and that have complex 
phases. We shall see that there is a corresponding nonlinear superposition 
principle for pure-state Wigner functions, that incorporates the relative phase 
of state vectors that are supersposed. 

2 Group representations 

Two examples for systems with one linear degree of freedom illustrate the 
strikingly different treatment of group representations in the Hilbert space 
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and phase space formulations [THj . 

Consider firstly the Heisenberg-Weyl group, which in the abstract may be 
considered as the 2-parameter abelian Lie group of transformations g(a,b), 
a, b e JR., acting on real, square-integrable functions F(q,p) of two real vari- 
ables as 

g{a, b) : F(q, p) — > F(q + a,p - b) , (2) 
We can write, on suitably smooth F, 

g(a, b) = expfaxi + bx 2 ] , g(a, b) g(c, d) = g(a + c,b + d) . 

x 1 = d/dq, x 2 = -d/dp, [x\ , x 2 ] = . (3) 

and we note the appearance of the 2-dimensional real Lie algebra spanned 
by X\ and x 2 . 

The familiar unitary ray representation of this group on Hilbert space has 
g(a,b) — ► U(a,b) , 

U(a, b)U{c, d) = e ^ ad - bc ^ 2h U{a + c,b + d). (4) 

With 

U(a, b) = exp[i(aP + bQ)/H] , (5) 

we have 

Q, P 

where I is the unit operator, so that at the level of the Lie algebra the ray 
representation of the group is associated with a central extension, charac- 
terised by Planck's constant h. In either the coordinate realization or the 
momentum realization of Hilbert space, one of the operators Q, P is real, 
the other pure imaginary, and again we see the essential role played by the 
complex numbers. 

Because this ray representation and central extension lie at the very heart 
of the Hilbert space formulation, it is all the more surprising to find that 
in the phase space formulation, the Heisenberg-Weyl group simply has its 
defining representation (j2JEJ). This representation is real and true; not only 
the complex numbers, but also Planck's constant, are nowhere to be seen. 
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(6) 



4 



The second example is provided by the time reversal group of order 2, 
with elements g, e acting on functions F(q, p) as 

g:F(q,p) — > F(q, -p) , e : F(q, p) — > F(q, p) , (=> g 2 = e) . (7) 

In the Hilbert space formulation, e is represented by the unit operator 
U e = I, and g by an antiunitary operator A g . In the coordinate realization, 
for example, we may have 

A g : ${x) — ► , Ue ■ il>{x) — ► V(ar) , (=S>Aj = l7 6 ). (8) 

Once again, the representation in the phase space formulation is simply the 
defining representation (J7j), and in sharp contrast to the Hilbert space repre- 
sentation, is real and unitary. 

It is not hard to find the source of these differences. The real structure 
that is inherent in the phase space formulation is in fact already present 
in the Hilbert space formulation. To see this, we recall jHj, ^E] that the 
mapping from the Hilbert space formulation to the phase space formulation 
is the Weyl-Wigner transform W, whose inverse W _1 is Weyl's quantization 
map. The transform W does not act on state vectors in Hilbert space at 
all. Rather it maps the real Hilbert space T of hermitian Hilbert-Schmidt 
operators A that act on the space of state vectors, onto the real Hilbert 
space /C of square integrable functions A on phase space. Conversely, W _1 
carries such functions A back into hermitian operators A (that is, it quantizes 
them). If Ak{x, y) is the kernel of A, regarded as an integral operator in the 
coordinate realization of Hilbert space, then the action of W and W" 1 is 
defined by 

W(A)(q, p) = J A K (q-y/2,q + y/2)e™/ h dy, 

W-\A) K (x,y) = A([x + y]/2,p)e^ x -y^ h dp. (9) 

With a suitably generalized interpretation of the integrals, these formulas 
also define the action of W on hermitian operators that lie outside T, such 
as Q, P, and the action of W -1 on real- valued phase space functions that lie 
outside /C, such as q and p. 

It can be seen from this that the question of the equivalence of the Hilbert 
space and phase space formulations of quantum mechanics is really a question 
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as to whether or not, in the usual Hilbert space formulation, one needs state 
vectors to truly describe states, or whether it is sufficient to work only with 
density operators and the real algebra of hermitian operators. 

The density operator p for the state of the system involved, lies in T and 
is mapped by W into (2Tchx) Wigner's pseudo probability density function 
W p G JC on phase space. But there is no mapping of the complex state vector, 
which has no image appearing in the phase space formulation. 

A ray representation it of a group G, defined by unitary and antiunitary 
operators acting on state vectors if), also defines a real, true representation 
n r of G on T, 

g GG : tjj — ► tt(#)^, 

A^U T (g)(A)=7r(g)An(gy. (10) 

It is this real, true representation Uq- that is mapped by the Weyl-Wigner 
transform onto a real, true representation il^ on square-integrable functions 
on phase space. 

From (jlUj) it can be seen that U/c is isomorphic to the tensor product of 
the ray representation ft with its contragredient ft c ; elsewhere [18 j we have 
called ILx the Weyl-Wigner product of fr and ft c . 

This explains why groups have real, true representations in the phase- 
space formulation. But if that formulation is indeed equivalent to the Hilbert 
space formulation, it must be possible to recover the underlying ray repre- 
sentations from the phase-space ones. In effect, we must be able to 'factorize' 
11^, or ITr, and determine ft. That this is possible is surprising, especially 
when it is recalled that at the level of the Lie algebra we have to determine a 
centrally extended representation, with associated parameter(s), from a true 
one. Because the representation LT^ can typically be regarded as the defining 
representation of the group of interest, as in the examples of the Heisenberg- 
Weyl group and time reversal group above, it can be said that in finding ft 
we are expressing Wigner's Theorem in a constructive way for that group. 

It has to be emphasized that not every real, true representation on K, of a 
given group can be factorized in this way. Necessary and sufficient conditions 
for this to be possible have been determined elsewhere JH] ■ We shall expand 
here on one aspect only: the construction of it in the case of a 1-parameter Lie 
group, whose representation U/c is generated by a real, skew-adjoint operator 
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on /C, which we may think of as an integral operator a with real kernel ajc, 
(&F) {q, p) = J a K {q, p, q', p') F{q', p') dq' dp' . (11) 
In this case, the necessary and sufficient condition on a for factorizability 



is 



sin[(gpi +pqi)/h] R(qi, Pi,q,p)dq 1 dpi 

sin[(gip + qip' + piq + p 1 q')/2h] R(q 1 , p h [q + q']/2, [p + p'}/2) dq% dp x 

- J sin[(gij9 - qip + p x q' - p 1 q)/2h] R(q u p u [q' - q]/2, \p' - p]/2) dq x dp x ,(12) 
where 

R(q, p, q', p') = a K {[q' - q}/2, [p' + p]/2, [q' + q]/2, [p'-p]/2) . (13) 

When this condition holds, we can obtain the (complex) hermitian generator 
A of the representation 7r on Hilbert space by setting 

A(q,p) = c + J sin[(qp + pq')/h] R(q' , p', q, p) dq' dp , (14) 

where c is an arbitrary constant, and then 

A = W-\A). (15) 

Note that A is defined only up to the addition of a multiple of / by c, leading 
to an arbitrary overall phase in the representation of the 1-parameter group 
on Hilbert space. 

It is surprising that complex ray representations can be obtained from 
true representations by this factorization process, which involves the 'quan- 
tization' step (|15p. At the level of the corresponding Lie algebra, we may 
obtain a representation of a central extension of that algebra in this way. 
Clearly any extension parameters must be built into the transforms (|T2"j) and 
W _1 (namely H), or else appear already in the true representation on phase 
space. Both possibilities arise in practice [TH] . 
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As an example, consider again the Heisenberg-Weyl group for one degree 
of freedom, with representation (0 EI) on phase space. In this case, as seen 
from (jSJ), cui = Xi — d/dq, a 2 = x 2 = —d/dp, leading from (fTTj) to 

<Xik(q, P, l', p') = S'(q-q^)6{p-p') , a 2K (q, P, q', p) = -S(q-q')5'{p-p') , 

(16) 

and hence from (|14j) and (|15|) to 

Ai(q, p) =p + ci, A 2 (q, p) = q + c 2 , 

A x = P + cj , A 2 = Q + c 2 I. (17) 

The constant terms in Ai, A 2 can be removed by a unitary transformation 
generated by Q followed by a unitary transformation generated by P, leaving 
P and Q satisfying (0). 



3 Superposition of Wigner functions 

According to the first of (jOJ), the Wigner function W p corresponding to 
([27r/l] _1 x) the pure-state density operator p with kernel px(x, y) = (p(x)<p*(y) 
is 

W p (q,p) = / <p(q-y/2)<p*(q + y/2)e^ h dy. (18) 



2-kK^ 

Here <f(x) is the representative of a normalized pure state vector correspond- 
ing to p, in the coordinate realization of Hilbert space. It is clear that <p(x) 
can be recovered from W p (q, p) by applying the second of Q- We have [TH] 

<p(x)<p*(x ) = [ W p ([x + x }/2, p) e^ x ~ x ^ h dp, (19) 



which determines a fixed complex multiple of <p(x) so long as we choose any 
particular x such that (p(xo) ^ and hence such that the left- hand- side of 
(I19jl is not identically equal to zero. What we wish to show here is that a 
complex superposition of two state vectors, including their relative phase, can 
also be described in terms of the two corresponding real Wigner functions. 

Suppose then that we are given two normalized, orthogonal pure-state 
wave functions, <fi(x), <f 2 (x), from which we have constructed the corre- 
sponding density operators pi, p 2 , and the corresponding Wigner functions 
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W%(q, p), W 2 (q, p). We consider the normalized superposition 



<p(x) = ai (fii(x) + a 2 <pz(x) , Wi\ 2 + |«2| 2 = 1 • 



(20) 



Choosing particular values X\, x 2 such that y?i(xi) 7^ 0, <p 2 (x 2 ) 7^ 0, we 
rewrite (f2T)|) as 



<f(x) = ci<fi(x)<fl(xx) + c 2 ip 2 (x)ip* 2 (x 2 ) , |ci| 2 |v3i(a;i)| 2 + \c 2 \ 2 \ip 2 (x 2 )\ 2 = 1 . 

(21) 



The kernel of the density operator corresponding to the pure state (p(x) is 
then 



Pk{x,v) = V>(x)ip*(y) = |ai|Vi(z)<pi(y) + \a 2 \ 2 tp 2 {x)ip* 2 {y) 

+c l c* 2 if l (x)ifl(x 1 )ifi* 2 (y)ip 2 (x 2 ) + c 2 cl(p- l (x)(pl(x2)ipl(y)ipi(x 1 ) , 

= K| 2 pi K (x,y) + \a 2 \ 2 p 2K (x,y) 

+c 1 c* 2 p 1K (x 1 x t )p 2K (y, x 2 ) + c 2 c\p 2K {x, x 2 )p* 1K (y, x x ) .(22) 
Applying the first of to p^/27r/l, we obtain 

W p (q, p) = laxl 2 ^^, p) + \a 2 \ 2 W 2 {q, p) 



Next rewriting pix and p 2 x using the second of (jUJ), we see that the third 
term on the RHS of ()23|) can be written as 






x e «Pi {q-y/2-xi]/K e -ip2 [q+y/2~x 2 ]/h e ipy/h ^ j ^ ^ 



(24) 
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where we have noted the reality of W 2 . Because the fourth term on the RHS 
of (|23|) is the complex conjugate of the third, we see that the sum of the third 
and fourth terms can be written as 



C l C 2 




W 1 {[q- y/2 + Xl }/2, Pl )W 2 ([q + y/2 + x 2 }/2, p 2 ) 



7rh 

x cos[e + (2p -pi- P2)y/2 + (pi - p 2 )q - p\X 1 + p 2 x 2 ] dpi dp 2 dy , (25) 
where we have written 

c lC ; = |ci||c 2 |e ie , 0<e<27r. (26) 

Next we note that 

M 2 = |ci|Vi(^i)| 2 = M 2 J W(x u p) dp, (27) 
with a similar relation for | a-2 1 2 - Then we can write (|23|) as 
W(q, p) = \c 1 \ 2 { J W^p^dp'jW^q, p) + |c 2 | 2 { J W(x 2 ,p')dp'}W 2 (q, p) 

+ \ciM rrr Wi ( [q _ y/2 + Xi]/2; Pl)w ^ q + y/2 + X2 y 2 , V2 ) 



irh 

x cos[e + (2p - pi - p 2 )y/2 + (p 1 - p 2 )q - p x x\ + p 2 x 2 \ dpi dp 2 dy .(28) 

This is a nonlinear superposition rule for Wigner functions, expressing 
a new Wigner function W in terms of two given ones, W\ and W 2 , two 
nonnegative coefficients |ci|, \c 2 \, and a phase angle e between and 2n as in 
(J2HJ). The two coordinates x±, x 2 in (|23j) are arbitrary, except that we must 
have 

y WiOr^dp^O, J W 2 (x 2 ,p)dp^0, (29) 

for consistency with our assumption above that <p(xi) ^ and <p(x 2 ) ^ 0. 
Since we also assumed that <fi and ip 2 are orthogonal, we must also require 
that Wi and W 2 are orthogonal, that is 

JJ W 1 (q,p)W 2 (q,p)dqdp=^\ J V \{x)ip 2 (x) dx\ 2 = , (30) 

10 



and since we assumed that W\ and W 2 correspond to pure states, they must 
also satisfy the pure-state conditions [TH] 

J W 1 (q,p) 2 dqdp=^, J W 2 (q,p) 2 dqdp=^ n , (31) 

as well as the normalization conditions 



/ 



W t (q, p)dqdp = l, J W 2 (q, p)dqdp = l. (32) 



It can now be seen that (j2*U|) and (|28|) are equivalent. If we are given two 
orthogonal, normalized wave functions tpx and tp 2 and form the superposition 
fOUJ) . then from ip we can construct W p as in (|28|) following the steps outlined 
above. Conversely, if we start with W\ and W 2 , we can find an X\ and x 2 such 
that (}2T?|) hold, and form the nonlinear superposition ()28|) for some choice of 
|ci|, \c 2 \ and e such that the second of (}26|) is satisfied. Note that fixing |ci|, 
| C2 1 and e fixes complex c\ and c 2 , up to an overall phase, as for example 
in (jnnj) below. We can also construct ipi(x)(p*(xi) and ip 2 (x)ip 2 (x 2 ) from W\ 
and W 2 , as in (JT5|) . and then recover ip as in (j2~Uj) . up to an overall phase. 

The result can also be expressed as follows. Given two orthogonal pure- 
state Wigner functions W\ and W 2 satisfying (jSOJ), (|3*T|) and (jHH), we choose 
any constants Xi and X2 such that (|29p hold, and any constants A, B and e 
satisfying 

A>0, 5>0, 0<e<2vr, (33) 

and we form 

T(q, p)=A 2 { j Wixup^dp'jW^q, p) + B 2 { J W(x 2 , p') dp'} W 2 (q } p) 
+ HlWxiiq- y/2 + Xl ]/2, Pl )W 2 ([q + y/2 + x 2 }/2, p 2 ) 



x cos[e + (2p -pi- p 2 )y/2 + (pi - p 2 )q - p x x x + p 2 x 2 \ dpi dp 2 dy .(34) 

Then 

W(q, p) = T(q, p)/ J T(q', p') dq 1 dp' (35) 

is a normalized pure state Wigner function. It corresponds to the wave func- 
tion superposition (}2~Tj) . where ipi(x)tp\(xi) and <p 2 {x)ip 2 (x 2 ) are constructed 
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from Wi and W 2 as in (|19p. and 

|ci|=A/{ y T(q, p)dqdp} 1/2 , \c 2 \=B/{ J T(q, p) dqdp} 1/2 , 

c\ = | ci | , c 2 = |c 2 | e~ le . (36) 

4 Conclusions 

Quantum symmetries and complex superpositions of pure, orthogonal quan- 
tum states, can both be described in the phase space formulation of quantum 
mechanics, without the use of complex numbers. The description of symme- 
tries is much simpler in phase space, but the description of superpositions is 
much more complicated. The description of non-orthogonal superpositions 
is even more complicated, and has not been attempted here. 

Note that the relative phase of the superposed wave functions, whose 
physical meaning is obscure [2|, is essentially the e appearing in the nonlinear 
superposition formulas (J2*%|) and (|34j). 

We conclude that the phase space formulation does indeed appear ca- 
pable of reproducing all aspects of quantum mechanics. In the case of the 
superposition of quantum states however, this is only be achieved at the cost 
of much greater complication. 

If we wish to think of the phase space formulation as the more fundamen- 
tal, arising directly from a deformation of classical mechanics in phase space 
|12| . we can think of the formulation of quantum mechanics in Hilbert space, 
and the associated introduction of complex numbers, as a computational 
device to make calculations easier. From this point of view, the apearance 
of complex numbers in quantum mechanics is on a similar footing to their 
appearance in other applications of mathematics to natural phenomena. 
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